Pants Decomposition of the Punctured Plane 
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Abstract 

A pants decomposition of an orientable surface E is a 
collection of simple cycles that partition E into pants, 
i.e., surfaces of genus zero with three boundary cy- 
cles. Given a set P of n points in the plane E 2 , 
we consider the problem of computing a pants de- 
composition of E = E 2 \ P of minimum total length. 
We give a polynomial-time approximation scheme us- 
ing Mitchell's guillotine rectilinear subdivisions. We 
give an 0(n 4 )-time algorithm to compute the short- 
est pants decomposition of E when the cycles are re- 
stricted to be axis-aligned boxes, and an 0(n 2 )-time 
algorithm when all the points lie on a line; both ex- 
act algorithms use dynamic programming with Yao's 
speedup. 

1 Introduction 

Surfaces (2-manifolds) , such as spheres, cylinders, 
tori, and more, are commonly encountered topolog- 
ical spaces in applications like computer graphics and 
geometric modeling. To understand the topology of 
the surface or to compute various properties of the 
surface, it is useful to decompose the surface into sim- 
ple parts. Among the possible ways to decompose a 
given surface, it is desirable to compute an optimum 
decomposition, one that minimizes a metric depend- 
ing on the application. 

A decomposition of an orientable surface E that 
has been studied is a pants decomposition E], a 
collection of disjoint cycles that partition E into pants, 
where a pant 1 is a surface of genus zero with three 
boundary cycles. Every compact orientable surface — 
except the sphere, disk, cylinder, and torus — admits 
a pants decomposition 

A natural measure of a pants decomposition to 
minimize is the total length of its boundary cycles. 
The length of a pants decomposition II of E, denoted 
by |IT|, is the sum of the (Euclidean) lengths of all 
the cycles in II. (If a subsegment is traversed more 
than once, its length is counted with multiplicity.) A 
non- crossing pants decomposition is a pants decom- 
position that allows any two cycles to touch as long- 
as they do not cross transversely. A shortest pants 
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decomposition is a non-crossing pants decomposition 
of minimum length. 

The problem of computing a shortest pants decom- 
position of an arbitrary surface E is open. In this 
paper, we study a variant of the problem where E is 
the punctured plane, i.e., E = E 2 \ P where P is a 
discrete set of n points. Figure ^i) gives an example 
pants decomposition of the plane with 6 punctures. 
Colin de Verdiere and Lazarus |2] studied a related 
problem: given a pants decomposition of an arbitrary 
surface, they compute a homotopic pants decompo- 
sition in which each cycle is a shortest cycle in its 
homotopy class. 





A A 



A A 

abed 



(i) 



(ii) 



Figure 1: (i) A pants decomposition of the plane punctured 
by a set of 6 points, and (ii) the corresponding binary tree. 

A cycle C on E is essential if it does not bound a 
disk or an annulus. A pants decomposition, also called 
a maximal cut system 0, is naturally obtained by 
the following greedy procedure. Let C be an essential 
cycle. Cut E along C to get a surface E' with two 
additional boundary cycles C\ and C% corresponding 
to the two sides of the cut. The cycles C\ and C2, 
together with the set of cycles obtained by recursing 
on E', is a pants decomposition II of E. The resulting 
cycle structure can be modeled as a binary tree with 
n leaves (Figure ^ii)). Each cycle C of II is essential 
because it encloses two other cycles, say C\ and C2; 
we write C\ -< C and C2 -< C to indicate that C 
encloses C\ and C2 . We say that two cycles C\ and C2 
are independent if neither C% -< Ci nor C2 -< C\. The 
final result is a pants decomposition of a bounded 
subset of the plane together with a single unbounded 
component. 



1 We call a surface of genus zero with three boundary com- 
ponents a pant instead of a pair of pants |2| . We refer to two 
pants instead of two pairs of pants; the latter phrase can be 
misunderstood to mean four such surfaces. 



In this paper, we give a simple algorithm with ap- 
proximation ratio O(logn), and a polynomial-time 
approximation scheme (PTAS) using Mitchell's guil- 
lotine rectilinear subdivisions. We compute the short- 
est pants decomposition in 0(n ) time when the cy- 
cles are restricted to be axis-aligned boxes, and in 
0(n 2 ) time when all the points lie on a line; both ex- 
act algorithms use dynamic programming with Yao's 
speedup, and are faster by a linear factor than the 
'naive' dynamic programming formulations. 



2 A simple approximation algorithm 

If II* is a shortest pants decomposition of the punc- 
tured plane £ = E 2 \ P, then every cycle in II is a 
simple polygon whose vertices belong to P. We argue 
next that II* contains a traveling salesperson (TSP) 
tour of the points. Choose any vertex on the outer- 
most cycle as the start of the tour. Traverse the outer- 
most cycle counterclockwise. The first time we visit a 
vertex u that also belongs to an inner cycle (one of the 
two legs of the pant) that has not been traversed yet, 
we recursively construct a tour beginning and ending 
at u that traverses the unvisited vertices on or in the 
interior of this inner cycle. After the recursion, we 
continue along the outermost cycle, repeating the re- 
cursive traversal of the second leg of the pant, until we 
reach our original starting point. Hence, II* must be 
at least as long as a shortest Euclidean TSP tour T* 
of the points. Hence, |n*| > \T*\. 

We convert a TSP tour T of the points in P to a 
pants decomposition n as follows. Initially, n is the 
empty set. Order the points from through n — 1 
in the order along the tour T. Let C(i,j) denote the 
polygon with vertices i, i + 1, i + 2, . . ., j — 2, j — 1, 
J: j ' ~ 1j J ' ~ 2, . . ., i + 2, i + 1, i in order, where 
the indices are taken modulo n. Imagine cycles of 
zero length around each point. Repeatedly introduce 
a new cycle into n' that is obtained by merging the 
two cycles adjacent along the tour enclosing the fewest 
number of points. Each cycle C(i,j) is obtained by 
merging two cycles C(i, k) and C(k + 1, j) by doubling 
the edge between vertices k and k + 1. We ensure 
that each edge in the tour T appears exactly twice in 
at most [logn] cycles in the pants decomposition H 
Hence, |H| < 2[logn]|T|. 

It is well-known 7 how to obtain a 3 /2-approximatc 
shortest Euclidean TSP of the point set using 
Christofides' algorithm in 0(n 3 ) time; the minimum 
spanning tree of the n points can be used to obtain 
a 2-approximation in 0(n log n) time. An approxi- 
mate TSP tour obtained by either of these algorithms 
gives us a non-crossing pants decomposition of length 
O(logn) times the optimum. 



3 PTAS 

Let e > be an arbitrary constant. To construct in 
polynomial time a (1 + e)-approximation to the short- 
est non-crossing pants decomposition, we modify the 
PTAS for Euclidean TSP tour due to Mitchell 
Our algorithm is more complicated because a pants 
decomposition consists of O(n) cycles instead of just 
one cycle as in a TSP tour. The PTAS is a dynamic 
programming algorithm where each subproblem is a 
rectangular region of the plane and two adjacent sub- 
problems interact only through O(l) grid points or 
portals. 

Let m > 2 be an integer and let M = m(m — 1). 
Let B denote the axis-aligned bounding box of the 
point set P. Imagine a shortest pants decomposi- 
tion n*. Mitchell [Hj has shown that there exists a 
favorable cut, i.e., a horizontal or vertical line I, which 
partitions B into two smaller boxes that can be recur- 
sively subdivided using favorable cuts. The recursion 
stops when a box is empty of points of P. Just like 
Mitchell, we introduce a segment of I, called a bridge, 
and O(M) grid points on I. Mitchell has shown that 
the total length of the additional subsegments is at 
most -^2 times the length of n* . 

Let R, a rectangle, be the boundary of an arbitrary 
box Q during the recursive subdivision. Intuitively, 
we can "bend" the cycles of n* to make each cycle 
that crosses R transversely do so only at one of the 
portals (grid points) and possibly use subsegments of 
the bridges on the four sides of R, without increasing 
the length of the pants decomposition by too much. 

To construct a short pants decomposition n, our 
PTAS solves subproblems of the following form. We 
are given a rectangle R whose sides are defined by 
two horizontal and two vertical favorable cuts. We are 
given two integers rii and n t , both in the range from 
through ri— 1, of the number of cycles of n that are in- 
side R and that intersect R transversely, respectively. 
Each of the rii cycles inside R intersects R tangentially 
and an even number of times, and each of the n t cy- 
cles intersects R transversely and an odd number of 
times. Let tir = rii + n t . We are given the pattern 
in which the tir cycles intersect R at the O(M) grid 
points on the sides of R. There are 0(n°( M )) possible 
ways for the tir cycles to intersect the sides of R. 

It remains to account for the fact that cycles in n 
that intersect R tangentially can traverse subsegments 
of the cuts bounding R. We observe that every point 
in the plane lies on at most two independent cycles 
of n. Let p be an arbitrary point in the plane. Let C p 
denote the subset of cycles in n that pass through p. 
If | C p | > 2, then there exist three cycles C\, C2, 
and C3 in C p such that both C\ and C2 are inside C3 . 
Let C be an outermost (minimum depth) cycle that 
traverses a subsegment ab of some cut. Each subseg- 
ment is shared by at most two independent cycles. 



Therefore, we count the length of ab at most twice 
when counting the total length of all subsegments of 
cuts traversed. 

The dynamic programming algorithm proceeds as 
follows. For a rectangle R intersected by tir = n,i + nt 
cycles, we try each of the 0(n) favorable cuts that par- 
tition R into two smaller rectangles, A and B. We try 
each of the 0(n^ M ^) possible ways that the tir cycles 
can intersect the cut transversely, making sure that 
the pattern in which the cycles intersects the cut is 
consistent with the pattern in which they intersect R. 
Some of the m cycles that belong inside R may belong 
inside A and some others inside B; we try the 0(ni) 
possible ways to allocate a subset of the rij cycles to A 
and the remaining to B. We optimize over the 0(n) 
cuts and 0(n°( M ^) intersection patterns to solve the 
subproblem R optimally. In the base case, if R has no 
points of P in its interior, then the subproblem has 
only 0(M) size, which is a constant, and is solved by 
brute force. Since there are 0(n°( M ^) different sub- 
problems and each subproblem takes 0(n°( M )) time, 
the total running time is 0(n°( M )). 

The length of the pants decomposition II ob- 
tained by the dynamic programming algorithm is 
o(l + times that of a shortest pants decom- 

position. To obtain the desired approximation factor, 
we choose m > 2\f2/e. 

To reiterate, the major differences between our 
PTAS and that of Mitchell are the following, 
(i) Each of our n — 1 cycles crosses transversely the 
boundary of a rectangular subproblem in one of O(M) 
different ways. Therefore, we have 0(n°( M ^) times as 
many subproblems to solve as in the TSP. (ii) Each 
of the 0(M) grid points on the boundary of a rectan- 
gular subproblem may lie on any of the n — 1 cycles 
in a pants decomposition. Therefore, the additional 
information associated with each subproblem is more 
than of constant size; the amount of information asso- 
ciated with each subproblem is 0(n°( M ^). However, 
the total running time is still polynomial in n. 

4 Points on a line 

Let P be a set of points on a line, without loss of 
generality on the £-axis. Order the n points from left 
to right. Let Xi denote the x-coordinate of the ith 
point. For every 1 < i < j < n, let (i,j) denote the 
consecutive points numbered from i through j. 
We prove next that a shortest non-crossing pants 
decomposition of £ = E 2 \ P must consist of convex 
cycles only. Hence, if n* is a shortest non-crossing 
pants decomposition of £, then there exists a A: in the 
range 1 < k < n such that n* consists of a shortest 
pants decomposition of (i, k) and a shortest pants de- 
composition of {k + 1, n) together with an outermost 
cycle of length 2{x n — X\) enclosing all n points. 



Suppose to the contrary that there is a cycle C 
in n* that contains a non-contiguous subset of points 
in P. Without loss of generality, we assume that C is a 
minimal one in the sense that both its legs, C\ and C2 , 
contains contiguous subsets of points. Without loss of 
generality, assume C\ is to the left of C2; thus, C\ is 
the left leg and C2 is the right leg of C. 

Let X C P be the set of all points between C\ 
and C2. Let x G X be arbitrary. Let T> x = 
{Di, D 2 , Ds, . . . , Di, . . .} be the set of cycles in n* 
containing x such that -Di+i contains Di. Let i be 
the smallest index such that Di contains some point 
of P \ X. There are two cases to consider: (i) Di 
contains C, (ii) Di does not contain C . 

If Di contains C, then we construct another cycle E 
enclosing C\ and D^_x making E the left leg of C 
(instead of C\), Delete Di. 

Otherwise, if Di does not contain C, then either 
(a) Di contains points only to the left of C2 or (b) Di 
contains points only to the right of C\. In the for- 
mer case, Di-\ is the right leg of Di. We swap C\ 
and Di-i so that C\ is the new right leg of Di 
and Di_i is the new left leg of C. In the latter case, 
Di-i is the left leg of Di. We swap C2 and so 
that C2 is the new left leg of Di and is the new 

right leg of C. 

In either case, we obtain a pants decomposition 
with total length smaller than n*, which is a con- 
tradiction. 

Let c(i, j) denote the cost of a shortest pants decom- 
position of (i, j). We have just proved that c(i,j) sat- 
isfies the following recurrence for every 1 < i < j < n: 

c (hj) = ~ Xi) + min (c(i,k) + c(k + 1, j)) (1) 

i<k<j 

where c(i,i) — 0. A shortest pants decomposition of 
(i,j) can be computed by choosing the appropriate 
value of k in the range i < k < j, computing the 
optimum pants decompositions of (i, k) and (fc+ l,j), 
and introducing a non-crossing cycle of length 2(xj — 
x^ enclosing the points (i,j)- The straightforward 
dynamic programming algorithm computes c(l,n) in 
0(n 3 ) time. 

We show how the running time of the dynamic pro- 
gramming algorithm can be improved by a linear fac- 
tor using Yao's speedup (HJ ■ Let w{i,j) = Xj— Xi. The 
function w() is monotone, i.e., w(i,j) < w(k,l) when- 
ever (i,j) C (k,l), and satisfies the following concave 
quadrangle inequality 

Vi < i' < j < j' : w(i,j)+w(i',j') < w(i',j) + w(i,j r ) 

In fact, the above equation is satisfied with equality 
because w(i,j) + w(i',j') = (xj — Xi) + (#'• — x^) = 
w(i', j)+w(i,j'). Let cu(i,j) denote w(i, j) +c(i, k) + 
c(k + l,j). Let K(i,j) denote the maximum k for 



which c(i,j) = Cfc(i,j). The following claims are al- 
most identical to those in the context of optimum bi- 
nary search trees proved by Mehlhorn |S] : 

1. The function c(i,j) also satisfies the concave 
quadrangle inequality, i.e., 

Vi<i'<j<f: c(i,j)+c(i',j') < c(i',j)+c(ij') 

2. K(i,j-l)<K(i,j)<K(i+l,j) 

We compute c(i,j) by diagonals, in order of increasing 
value of j — i. For each fixed difference d, we compute 
c(i,j) where j = i + d; we compute Ck(i,j) for k in 
the range K(i,j — 1) < k < K(i + 1, j). The cost of 
computing all entries on the dth diagonal is 

n — d 

J]jf(i + i,j)-jr(tj-i) + i 

= K(n - d + l,n+ 1) - K(l,d) + n-d 

< (n+l)-l + n-d 

< 2n 

Since d ranges from through n — 1, the total running 
time is 0{n 2 ). 

5 Box decomposition 

A box decomposition of E is a pants decomposition II 
in which each cycle in IT is an axis-aligned rectangle. 
Observe that any two axis-aligned rectangles can be 
separated from each other by either a horizontal or a 
vertical line. 

Let X\ through x n denote the ir-coordinates of the n 
points in increasing order, and let y\ through y n de- 
note the y-coordinates of the n points in increasing 
order. Let h(i, j) = 2(xj — Xi) and let v(i,j) = 

2 {yj~Vi)- Let w(ii,i 2 , ji, 32) = h(ii,i 2 ) + v(ji,j 2 ); 
then, w(i\, i 2 , ji, 32) is the perimeter of the axis- 
aligned box whose sides have x-coordinates x% t and 
Xi 2 and y-coordinates yj 1 and yj 2 . The function w() 
is monotone and satisfies the concave quadrangle in- 
equality. 

The rest of the proof is similar to the case for 
points on a line. Let c(ii, i 2 , ji, ji) denote the cost 
of a shortest pants decomposition of the points in 
the box [x n ,x l2 ] x [y n ,y n ]. The cost c(h, i 2 , ji, j 2 ) 
obeys a recurrence that is a two-dimensional general- 
ization of Equation ^ Similar to the dynamic pro- 
gramming algorithm for points on a line, we com- 
pute c(i\, i 2 ,ji,j 2 ) by diagonals, in order of increasing 
value of max{i2 — ii, j 2 — j — 1}. For each pair of dif- 
ferences d\ and d 2 , we compute 0(^1,12,^1,^2) where 
ii = h + di and 32 = ji + e?2- The cost of comput- 
ing all entries on the diagonals defined by (d±,d 2 ) is 
0(n 2 ); since there are 0(n 2 ) such pairs (dx,d 2 ), the 
total running time is 0(n 4 ). 



6 Work in progress 

We mention some very interesting open questions that 
we are currently investigating. 

Is it NP-hard to determine, for an arbitrary L, 
whether there exists a non-crossing pants decompo- 
sition of the punctured plane of length at most LI 
Is there a simple algorithm to compute an 0(1)- 
approximate shortest pants decomposition? 

Are the cycles in a shortest (non-crossing) pants de- 
composition always convex? If not, how much longer 
than optimum is a convex pants decomposition? 

How efficiently can we compute a shortest pants de- 
composition of the plane with different types of punc- 
tures, e.g., rectangular holes instead of points? 

How efficiently can we compute a shortest pants 
decomposition of other 2-manifolds, such as the torus 
minus a set of points? 
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